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Abstract 



We give a proof of transient fluctuation relations for the entropy production (dissipation function) 
in nonequilibrium systems, which is valid for most time reversible dynamics. We then consider the 
conditions under which a transient fluctuation relation yields a steady state fluctuation relation for 
driven nonequilibrium systems whose transients relax, producing a unique nonequilibrium steady state. 
Although the necessary and sufficient conditions for the production of a unique nonequilibrium steady 
state are unknown, if such a steady state exists, the generation of the steady state fluctuation rela- 
tion from the transient relation is shown to be very general. It is essentially a consequence of time 
reversibility and of a form of decay of correlations in the dissipation, which is needed also for, e.g., 
the existence of transport coefficients. Because of this generality the resulting steady state fluctuation 
relation has the same degree of robustness as do equilibrium thermodynamic equalities. The steady 
state fluctuation relation for the dissipation stands in contrast with the one for the phase space com- 
pression factor, whose convergence is problematic, for systems close to equilibrium. We examine some 
model dynamics that have been considered previously, and show how they are described in the context 
of this work. 

Keywords: nonequilibrium phenomena, entropy production, chaotic hypothesis, time reversibility, 
correlation decay 
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1 INTRODUCTION 



Steady state fluctuation relations provide information about the probability distribution of time aver- 
ages of certain phase variables for systems in nonequilibrium steady states: they express the ratio of 
the probabilities of observing positive and negative values of time averages of such variables. These 
relationships have received much attention over the past decade (see e.g. [TJ EJ [3J Hj ) . In fact, they are 
among the few exact results regarding nonequilibrium systems; close to equilibrium they lead to the 
Green-Kubo relations for linear transport coefficients [5]; they provide information about fluctuations 
in nanoscale systems and they show how irreversibility can emerge from time reversible dynamics. 

The first such fluctuation relation (FR) for a nonequilibrium particle system was given in Ref.[6j, 
in 1993. It concerned fluctuations of time averages of the energy dissipation or equivalently the phase 
space contraction rate, for an isoenergetic shearing interacting particle system. In the isoenergetic 
dissipative system of Ref.[B], the instantaneous energy dissipation equals the instantaneous phase space 
contraction rate. The form of the relation was the following: 



Pr(A) 
PA-A) 



where A and —A were the time averaged values of the dissipation on steady state trajectory segments 
of duration r, and P T (A) is the steady state probability of observing the average value A (with some 
tolerance), in a time r. In analogy with the periodic orbit expansion for dynamical systems [7], 
the relation was derived heuristically using the "Lyapunov weights" in the long r limit, and was not 
expected to work at short r's. 

Paper [H] motivated a number of works, in which various derivations of formulae formally similar 
to Eq.fJIJ were given for different systems, although for non-isoenergetic systems the phase space 
contraction rate and the energy dissipation are no longer instantaneously equal. Beginning in 1994, 
a series of papers by Evans and Searles focussed on the fluctuation properties of the "dissipation 
function" D,, cf. Refs.pQ El El HD1 El H21 H31 E] , while Gallavotti and Cohen, in 1995 et.seq., considered 
the fluctuations of the time-averaged phase space contraction rate, —A, under the assumption that the 
dynamics is time reversible, transitive and Anosov or that the system satisfies the Chaotic Hypothesis 
(i.e. although the dynamics is not Anosov, for the purposes of this FR, it behaves as if it was Anosov), 
cf. Refs.[H3 [TBI HZ]- Various generalizations and extensions of these relations have been produced by 
different authors, see e.g. Refs.[TOJ EH 03 EDI ED EZl E3] • 

In the literature, FR's which superficially resemble the form of Eq.Q, are often collectively referred 
to as the Gallavotti-Cohen Fluctuation Relations, or instances of the Gallavotti-Cohen symmetry, 
although they may have subtly different meanings. Therefore, for the sake of clarity, we refer to the 
relation for the phase space contraction rate of Refs.jEJ Q5], as the A-FR, while we refer to the relation 
for the dissipation function of Refs.pQ EJ [5J H], as the f2-FR. In the cases in which = —A, like those 
of [§], the two relations have the same content. 

J7-FR's have been widely used and tested numerically and experimentally, for various transient 
[HIE] as well as steady state systems (see e.g. [3 [HI [lU HH1 E31 EH 123 EHl ) • Although the transient 
O-FR's have been proven more than a decade ago, beginning with Ref.jEJ E], there has been some 
question over the derivation of the corresponding steady state relations [25], Therefore, we present 
a derivation of the O-FR's which is more thorough than the ones that have appeared thus far in the 
literature [JJ. The purpose is to shed light on the mechanisms which allow the O-FR's to hold, and 
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which may be used to better understand properties of nonequilibrium systems. In particular, we are 
interested in nonequilibrium steady states where the dissipation function takes the form 

Q = £ = -JVF e /{k B T) , (2) 

where J is the dissipative flux in the system of interest, V is the volume of the system of interest, 
ks is Boltzmann's constant, T is the equilibrium thermodynamic temperature of the thermal reservoir 
with which the system of interest is in contact, and £ is the generalized entropy production (or energy 
dissipation) induced by the dissipative field F e applied to the system of interest [HE]. For the equality 
= £ to hold, we require that the unthermostatted (i.e. adiabatic) field dependent dynamics preserves 
phase space volumes instantaneously, and that the initial distribution of phases could in principle be 
generated by a single (exceedingly long) zero-field phase space trajectory. Systems that satisfy the 
condition of the adiabatic incompressibility of phase space are called AIT in Ref.(2Sj (cf. Appendix A). 

When referring to steady state systems, we only consider systems which relax to a unique steady 
state. The necessary and sufficient conditions for the system of interest to relax to a unique steady 
state are unknown. Certainly, it is known that under certain circumstances nonequilibrium steady 
states are not possible [23] • If the systems of interest are near equilibrium, where the thermodynamic 
notions of temperature and entropy can be used, then as is shown in the theory of linear irreversible 
thermodynamics the time or ensemble average of fl is equal to leading order in the dissipative field, to 
the average so-called spontaneous entropy production rate [31] • Outside this local equilibrium regime, 
the entropy production rate cannot be defined. While the notions of a thermodynamic temperature, 
thermodynamic entropy or an entropy production rate cannot be defined outside the local thermody- 
namic equilibrium regime, in a nonequilibrium steady state the ensemble or time averaged phase space 
contraction rate is equal to the steady state time derivative of the fine grained Gibbs entropy. The 
fact that the fine grained Gibbs entropy is not constant in a steady state (unlike the pressure or the 
energy), is the reason why it is so difficult to attach meaning to standard thermodynamic quantities 
such as temperature, in steady states outside the local equilibrium regime. 

In the derivation of a relation like ([!]), one may follow different approaches. The one which has 
led to the A-FR requires the full knowledge of the physical measure of the dynamics, which is then 
assumed to be of the Anosov type, hence to have a Sinai-Ruelle-Bowen (SRB) measure, hsrb [13] 152]. 
This is the mathematical approach, which aims to identify the class of dynamical systems for which 
a relation like Q can be rigorously obtained, regardless of the physical relevance of such dynamical 
systems or of the quantities appearing in Q. How far one can go along this line and which other 
condition could replace the the Anosov condition are highly nontrivial and interesting mathematical 
questions, which do not have an answer, at present. 

The explicit knowledge of the physical measure gives the maximum possible amount of information 
about the phase space distribution: all statistical properties of the dynamics can be obtained from 
it. However, if one is only interested in why the FR holds, less detailed knowledge of the phase space 
distribution is likely to suffice, and special features of the dynamics, such as the Anosov property, 
need not play a role. This may explain why many physical systems have an observable which obeys a 
relation like 0. As a matter of fact, hardly any particle system of physical interest is of the Anosov 
kind, therefore it is interesting to ask whether a physical system obeying a fluctuation relation does so 
because it shares certain properties with Anosov systems, or whether it obeys such relations for other 
reasons. To investigate this question, which is part of the physical approach, we follow Evans and 
Searles [IT], since they do not require the invariant phase space probability measure to be explicitly 
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determined. This will lead to transient fluctuation relations for S as well as for many other functions 
of phase, including A. If cj) is the observable of interest, we speak of </>-FR. 

The purpose of the present paper is to understand why the S-FR holds for many physically inter- 
esting systems, and why this can be observed within physically relevant time scales. This is useful to 
further develop the present theory. 

Remark. This issue is not completely unrelated to, but is quite different from the identification of the 
class of dynamical systems which allow a rigorous mathematical derivation of the A-FR. One may say 
that the two issues are as much related as Khinchin's approach to the ergodic hypothesis pJity is related 
to ergodic theory. 

In Section 2 we derive transient FR's for time reversal invariant dynamics. These relations are called 
"transient" because they describe the time evolving statistics of an ensemble of finite time averages of 
given phase functions. They hold under very general conditions, similarly to other transient relations, 
like the Jarzynski and the Crooks relations [34J. In particular, the transient relations hold even if 
a steady state with positive and negative fluctuations of the observable under consideration is not 
reached. If a fluctuating steady state is not reached, the transient relations hold and express properties 
of evolving ensembles of phase space trajectories, while the Anosov structure cannot characterize the 
dynamics. 

In Section 3, we show how the results of Section 2 can be used to derive asymptotic FR's for 
reversible systems, under the technical assumption that O is bounded. In Section 4, the decay of 
the ^-autocorrelation is shown to be sufficient for the steady state O-FR to hold. Furthermore, in 
isoenergetic systems which enjoy AIT (e.g. systems whose adiabatic evolution is Hamiltonian), the 
phase space contraction rate and the dissipation function are instantaneously equal. In these cases, 
the predictions of the A-FR and of the fi-FR can be compared directly, and the Q-FR at high dissipation 
implies a different behavior from the one predicted by the A-FR under the axiom-C condition of Ref.|35J, 
but which agrees with the results of [HI This is possible because the conjectures of [33 ES] are 
not necessarily verified by particle systems. The approach followed in Sections 2, 3 and 4 leads to a 
number of 0-FR's and to various other relations. In Section 5 we discuss our results. Appendix A 
contains a detailed description of the dissipation function for the common cases. Appendix B discusses 
the problem of axiom-C systems and the decay of correlations, and explains that our approach is based 
on exact relations, which are then valid even when the conditions of Section 4 are not met. Our 
conclusions are that: 

• to satisfy the steady state fi-FR, the time reversal invariance of the dynamics and the decay 
of the ^-autocorrelations, with respect to the initial (non-singular) probability measure, suffice. 
This decay of correlations is not the one characterizing the steady states of axiom-C systems, 
which concerns all observables and is referred to the invariant (singular) measure. 

• For f2 to equal E, hence to obtain a steady state S-FR, AIT must hold and the initial distribution 
of phases must in principle be obtainable form a single field free (F e = 0) phase space trajectory. 
In such a case, we say that the initial distribution is ergodically consistent with the field free 
dynamics [I]. 

• The fi-FR's avoid the difficulties encountered by the A-FR close to equilibrium. The A-FR, 
requires ever longer convergence times, the closer the steady state is to equilibrium (TTJ EU ESI 
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157] . The O-FR's hold in a range which does not shrink for decreasing fields, and are verified within 
the physically relevant (^-autocorrelation) times, both for systems far and close to equilibrium. 

• The direct derivation of the steady state S-FR's avoids the difficulties connected with obtaining 
a relation for the energy dissipation via the A-FR, justifies the convergence rates, and leads to 
new testable relations. However, the class of dynamical systems which enjoy the required decay 
of correlations remains to be identified. 

• It goes without saying that the steady state relations require that a steady state is reached (hence 
that the correlations decay with respect to the initial state), and that fluctuations of opposite sign 
in the dissipation function should be observable. If such a steady state is not attained or if the 
steady state is sufficiently far from equilibrium that negative values of the dissipation function 
cannot be observed, the steady state FR's cannot be applied, but the transient FR's (including 
their asymptotic forms) remain valid as properties of the evolving ensembles, because they only 
require time reversibility. 

In this sense, our relations are like the usual thermodynamic relations which are extraordinarily general, 
being independent of the nature of interparticle interactions and dynamics. All that is required of the 
dynamics is that it should be time reversible and, for steady state FR's, that the autocorrelation of 
O decays and that the steady state is uniqu^j] Furthermore our derivation explains why the O-FR is 
verified within times which can be expressed in terms of material properties. 



Because we will be interested in thermodynamic systems which possibly relax to a nonequilibrium 
steady state, our dynamics will consist of the Hamiltonian dynamics of an A^-particle system of in- 
teracting particles to which a dissipative field is applied. The resultant equations of motion may or 
may not be derivable from a Hamiltonian [55] • In order for such a system to be capable of relaxing 
to a nonequilibrium steady state, the system must be allowed on average, to loose heat. This can be 
accomplished by surrounding the system of interest with thermostatting walls which serve two pur- 
poses: to confine the system of interest and also to remove dissipative heat generated (on average) by 
the dissipative field applied to it. It has been shown that any loss of heat from a Hamiltonian system 
requires that the dynamics over the degrees of freedom in the system of interest is no longer volume 
preserving [3Sj. These details are explained further in the Appendix. In what follows however, we 
will consider a more general setting that is time reversible and does not necessarily preserve the phase 
space volumes. 

Given the phase space M of a particle system, consider a probability measure juonX, with density 
/, i.e. d/i(r) = f(T)dT for every point T G M. The measure \i does not need to be produced by any 
dynamics on Ai, although it could represent, for instance, the initial equilibrium distribution of the 
system. Further, given a (sufficiently well behaved) phase function : M. — > M, the probability that 
it takes values in the interval (a, b) C M, according to /i, is given by 



2 TRANSIENT FLUCTUATION RELATION 




(3) 



lr rhis is analogous to the equilibrium thermodynamic requirement of a unique allotrope. 
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where 

^■^{T E M : cf>(T) e (a,b)} , (4) 

is the set of points of the phase space for which (ft takes values in (a, b) J^] 

Let S T : M. — ► M. be the time evolution operator, which takes any point T G to its evolved 
image S T T, under some dynamics applied for the time r J^] We refer to time reversal invariant dynamics, 
i.e. dynamics which obey 

iS T T = S~ T iT for all T € M and all r G M (5) 

where i : M — ► .M, which obeys ii = i 2 = identity, is an involution representing the time inversion 
operator (e.g. ir = i(q, p) = (q, — p)) for the dynamics. Consider the phase functions defined by: 

Kt 0+ r(n=- f ° +T HS s T)ds , & 0lt0+ r(r):=r^ 0it0+r (r) , for t ,TGiR; reM (6) 
r -/to 

Take <5 > 0, to = 0, introduce the sets Aj = (A - 5, A + 5), Aj = (—A - 8, -A + 5) and, for an odd 
phase variable (i.e. 4>(iT) = — (/>(T)), consider the ratio 

i.e. the probability according to /i that 4>q t & A"g, divided by the probability according to the same \i 
that c/> G Aj. 

To compute this quantity, observe that the points of M. which fall in 4> 0t \ a - are those, and only 
those, obtained by doing the time inversion of the evolution, for a time r, of the points in O T L + , i.e. 

' S 

<Ao,tIa7 = iST h,r\ A + ■ ( 8 ) 

o o 

Indeed, take any T G <j) 0r \ A -, invert it and evolve it backward for a time r; the resulting point 
Y = S~ T iT is the initial condition of a trajectory segment of duration r over which one obtains 
4>q^ t G Ag, hence Y G 0o,rU +- Moreover, there are no points of iS ,r ^> 0iT | j 4+ which do not lie in </>o,tIa~ ■ 
This allows us to compute the denominator of Eq.(7| in terms of 4> 0t \ a +, through the coordinate 

1 5 

transformation 

dT 



T = iS T X , whose Jacobian is J 



dX 



exp ^J\{S s X)ds^j = e Ao ^ x) . (9) 



Here, the quantity A is the phase space expansion rate (the opposite of the contraction rate) which, 
for dynamics T = G(T) on M, is defined by 

A(r) = v-r = v-G(r) . (10) 

2 4>\(a,b) is assumed to be /i-measurable, like all sets in this paper. 

3 In particular, the dynamics could be determined by the equations of motion of a particle system subjected to a 
dissipative field F e , and to a deterministic thermostat meant to remove the excess energy pumped into the system by 
the field, so that a nonequilibrium steady state can be reached [29] ■ In this case, V — (q,p) denotes the generalized 
coordinates and momenta of all the particles comprising the system. 
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This leads to: 



/ f(T)dT= f f(iS T X) ^- dX = I f(iS T X)e A ^(X) d X (11) 

J *0,r\ A - J ^,r\ A + dX J <t>0,r\ A + 

5 6 A S 

A necessary condition for the derivation of the transient fluctuation relations is that 
/ f(T)dT + for all A for which / f(X)dT + 

6 5 

This condition is satisfied by appropriate choices of fi, for the given S T . For example, [i could be 
selected to be the probability distribution generated by the field-free (equilibrium) dynamics of a 
particle system, Sq say, for which S T is the nonequilibrium dynamics; i.e. \x could be generated by 
Sq , obtained by setting the dissipative field of the nonequilibrium dynamics S T to zero. In fact, this 
selection ensures that f(iS T T) / whenever f(F) + oQ Then, assuming that / is also time reversal 
invariants (i.e. that f(T) = f(iT) for every T G M.) the ratio of probabilities (uj becomes: 



= « = i 12) 

M<Ma 5 ") k , T \ A+ f(iS T X)exp(A 0;T (X))dX J^ L+ f(S-X) exp (A , T (X)) dX 



The restriction on the choice of \x for the given S T , i.e. that f(iS T T) / whenever f(T) / 0, is referred 
to as ergodic consistency of / with S T [I]. The main quantity used below can now be introduced. 

Definition. The time averaged dissipation function flt ,t +T '■ «M — ► M, for a time reversal invariant 
phase space probability density f is defined by /T^: 



i rto+r 

n t0 . t0+T (ry.=- / n(S s T)ds 

T J t 



provided f is ergodically consistent with S t . The instantaneous dissipation function is obtained 

from differentiation with respect to r of&o tT , for sufficiently regular f. 

Here the condition of ergodic consistency plays an important role, although it may not be so evident 
in general. Indeed, given a particle system with certain conserved quantities under its equilibrium 
dynamics Sq , the nonequilibrium dynamics S T may not conserve the same quantities, and a trajectory 
S T T may wander outside the support of /i, making the corresponding Q unusable in the following 
derivations. Thus, in physical applications / must be chosen with care, although many phase space 

4 Here, it is assumed that equilibrium dynamics corresponds to no net dissipation, and that typical trajectories explore 
all the phase space. 

5 The assumption that / be time reversal invariant is not necessary for our main results. A different choice will simply 
result in a alternative definition of the dissipation function, [T]. However, if /i is selected to be an equilibrium measure, 
then this condition is guaranteed, and relying on it makes the calculations more elegant. 



1 

T 
1 

r 



In 
In 



/(■s*°r) 

fiiS^Y) 
f(S t0 T) 

/(5*o+^r) " 



to+T 



A(S s T)ds 



'to 

^■to,to+r(r) 



(13) 
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densities are acceptable from a mathematical point of view. Using Eq.(13), Eq.(12| implies what has 
been referred to as the 



Transient cft-FR: 



M<?WL+) Uoj^ 

/~7 I ~ \ = 7 r o Tyv flirty := < exp (~ n ^)h eA+ (14) 



Here (exp (— Qo^))^ &A + is the ensemble average of exp[— Qo tT (X)], over the set of trajectories which 

satisfy the constraint that </> 0jT G ^ [TJ [TU]. In the case that (ft 0T = Oo, r the fluctuation relation 
assumes the particularly elegant form, called 



Transient ft-FR: 



M^q,tU+) _ _ [A+£(MjT)]r 



Here, e is an error term of magnitude \e\ < 5, which depends on 5, A and r, and appears because f2o,r 
is not necessarily constant and equal to A. Therefore, e(5, A, r) — > as 5 —* 0, for all ^4 for which the 



ratio of Eq.([T5]) exists [BE1E]. 

Remark. Different choices of f are possible, which lead to different £l-FR 's. In particular, equilibrium 
probability densities of many-particle systems which obey AIT lead to an Q-FR which concerns the 
entropy production, or the energy dissipation rate £ (cf. Appendix A). The uniform density in a compact 
phase space, f(T) = say, yields O = —A. A density such that fdT approximates as accurately 

as desired the (possibly singular) steady state measure is also allowed. 

The above relations are called "transient" because, independently of the length of r, they express 
properties of the measure /i and not of the possible steady state. The time r, indeed, enters the </>-FR 
only in the definition of the observable (ft T , and not of the phase space probability density. Equations 



( 14|15 ) are exact and do not need any limit to be taken either on r or on 5. As the calculation shows, 



the only properties that are required for Eqs.( 14|15 ) to hold are that the dynamics S T be time reversal 



invariant, to ensure that, for odd (ft, —A can be observed if A can, and that f(iS T T) ^ for any T for 
which /(r) / 0. The range of A's that can be observed depends both on S T and (ft. 

If the initial distribution can be generated by a single field free phase space trajectory, and AIT 
holds, the dissipation function O takes the form given in ^ and at equilibrium (where F e = 0), = 0, 
hence 

^Ml = (e oyJ + = 1 (16) 

o 

as it should be for the equilibrium measure fi, and for any odd variable (ft (for instance, (ft could be the 
instantaneous current J, which does not need to vanish at equilibrium, although its mean certainly 
does). For symmetric intervals (—5,6), one also has A = and can write 

ufaA-M) = ! hence (e -n , T) l (17) 
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In the limit that 5 tends to infinity, the second equality in (17 1 tends to a full phase space average, and 
produces the NonEquilibrium Partition Identity (also referred to as the Kawasaki normalization) [5| |2"9"] . 
This identity can be used to test the accuracy of the numerical simulations, and has even been used 
to calibrate experimental equipment [26J. The highly asymmetric convergence of numerical estimates 
to unity means that estimating the statistical uncertainty of experiments is problematic without this 
relationship [5]. If the dissipation function is bounded (i.e. |0| < fi*, for some f2* > 0), one obtains 

e -m- < ggvkj < e «r , (18) 

o 

for all odd 4>- Because different dissipation functions are allowed and the above ratio holds for all 
of them, a sharp bound can be given for all (p and independently of the physically relevant dissipa- 
tion function, by taking the smallest possible Q*. This constitutes a new prediction, which may be 
interesting to test in nonequilibrium systems. 



3 ASYMPTOTIC FLUCTUATION RELATIONS 

In order to derive steady state fluctuation relations, one may develop further the transient </>-FR's 
obtained in Section 2, considering the time averaged phase variable 4> t(hto+T (cf. Eq.(|6j)), where the 
time averaging begins at a time to > 0, rather than at time 0: 

5 

Here, 4> t0tto+T \E = {T e M : 4> tQttQ+T (T) G E} is the set of all phase points which are initial conditions 
of trajectory segments of duration to + r, over which 4> tQ j Q+T takes values in E, and E = Af or Aj . 
The calculation of the ratio (19) requires that the points of M lying in 4> to to+T \ A - be identified. To 
do that, consider a trajectory of length t = 2to + r, and take W in </> to to+T | ,+. Then, let X = S to W, 

Y = S to+T W, Z = S 2to+T W = S l W, and V = iZ = iS l W (see Figure 0. The trajectory segment 
between iY and iX produces the opposite of the average produced between X and Y, in the time 
interval (to, to + r), so that 

^to.to+rU- = iS % ,t +r\A+ ( 2 °) 

and, as in the previous section, 

/(T)dT = / f(iS l W) "'' 



r 

jd> 



t 0. t 0+ rl A 



dW 



dW 



[ /(S^e^iAoAW^dW (21) 





where J = \ -^ \ = exp[Ao,t(VF)] and we use time reversal invariance of /. 
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Figure 1: Schematic diagram illustrating the time-reversible evolution of two points in phase space 
that are related by a time reversal mapping, i. If the average of <\> between X and Y is A, the average 
between iY and iX is —A. 
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Note that the integral defining J is along the full trajectory of duration t, although cp to tQ+T 
only one of its parts of duration r. From the definition of Qot (Eq.(13)), one obtains: 



concerns 



A ,t(W) = In ~ Oo,t(W) 



f{S l W) 



and 



*0> t 0+ T l 



f{T)dT 



<t>t ,t Q + T \ A + 



f(W)e- n °^ w) dW 



(22) 
(23) 



The probability ratio of Eq.(19) can then be written as: 



K<f>; 



to,to+T\Af' 



t +T\AZ 



h i f(W)e- Qo ' t( - w Uw' 

J \,to+r\ A + J V > 



Id,. 



f(W)dW 



-i 



(exp (-fi 0) t)); 



"t ,t 



(24) 



which is the inverse of the average of exp( — Qo,t), conditioned to the constraint cj>f ^ _^ T € Ap, and the 
special case 4> = Q, yields 



to+T\A 



' ( -(exp(-Oo,t)) nt \ + eA 



O 

The right hand side of this expression can be written as: 



(exp (— f2 ,t)>7 



■to>*o+ 



e [A+£(5,t () ,A,r)]r / e -n ,i -nt + r ,2i +T\_ 1 

\ /^i ,i + T 



6A^ 



(25) 



(26) 



Taking the logarithm, and dividing by r one then obtains: 

A + e(5,t ,A,T) - iln( e - n °.*o- n 'o+r,2t 0+ r) 



1 M(^to,to+-rU+] 

-In— = s — 



(27) 



where the correction term e depends on 5, to, A and r and, like in Section 2, |e| < 5. This result is 
exact, and holds for all to, t and (5, and for any A for which both the probabilities in the left hand side 
of (27) do not vanish. It rests only on the time reversibility of S 1 *, and on the ergodic consistency of / 
with SK 



Equations (24) and (27) are valid for all r, hence also in the limit r — » oo, but remain transient 
relations because no matter how large r is; they refer to the initial measure (j,, even if the dynamics 
makes \i evolve into \i T . However, the above derivation suggests that for systems that reach a steady 
state, it may be possible to derive steady state FR's from the transient FR's. The idea is that for 
sufficiently large times to, Eqs.( 24|2"7 ) should concern trajectory segments of length r distributed 
according to the invariant measure, or very close to that. Furthermore, taking r long, and 5 small, 
should make the last two terms in Eq.(27) negligible. 



To investigate this possibility, we firstly transform Eqs.(24) and (27) by transferring the time 
evolution up to to from the sets of phase space points (f> to to+T \ ( m , to the phase space probability 
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distributions. Given an initial phase space probability density /o, and the corresponding probability 
measure fiQ, we denote their evolved counterparts after a time to by ft , and //t respectively: 



m {S t °E) = ii Q {E) or equivalently / f to (X)dX = [ fo(W)dW 

JS^E JE 



(2f 



for every to € M and E C M, where S to E is the set of points reached by all points of E after 



an evolution of duration to- Equations (28 1 are consequences of the conservation of probability by 
time evolution]^] Physically, they are an instance of what is referred to as the Heisenberg-Schrddinger 
equivalence, while mathematically they are the result of the coordinate transformations X = S to W 
and W = S~ to X. The Jacobian of this transformation is given by 



dW 



dX 



3 A ,_ t0 (X) = L *° A(S°X)ds 



-A_t ,o(X) 



(29) 



hence one can write 



t Q ,t +T\(a,b) 



and 



^o(^t ,t +r\{a,b)) 

b)) = ^t {S t0 (/) t0tto+T \ {afi) ) : 



fo(W)dW = 

f {S- to X)e- k - t O' o{x) dX 



(30) 



ft (X)dX . (31) 

'S*°<^ ,t Q -Hrl(o.b) 

In Eqs.( p0pl| ), we could have integrated over any measurable set E, hence, equating ( |30| ) and ( |3l| ) 
finally gives: 

f to (X) = f (S- to X)e- A -*o,o(*) . (32) 
This is known as the Lagrangian form of the Kawasaki distribution function [S] Q From Eq. ( 22 ) , which 

-to / f{X ) 



yields 



A- to ,o(X) = [ °MS s X)ds 
Jo 



In 



/(s-*°x; 



n 



,-t 



(X) 



Eq.(32) can be rewritten as 



f t0 (X) = /o(I)e- n ».-'oW 



(33) 



(34) 



and gives the evolution of the phase space probability density fo- This well known exact result for 
the time dependent phase space density ft (X), was first derived for thermostatted particle dynamics 
where the dissipation function takes the form given in Eq.(|2]) by Evans and Morriss in 1984 |39J. 
Observe that 

S to 4> t0)to+T \ {afi) =0 o ,tI(o,6) ( 35 ) 



6 They merely state that the initial probability of E becomes the probability of the evolved set S t<s E after a time to. 
7 This can also be written, ft (S to W) = f (W)e- A °- t o < - w \ 



13 



because 4> is a phase variable, hence does not explicitly depend on time, and because we assume the 
dynamics to be autonomous. Indeed, a point X S S to 4> t0tto+T \^ a ^) is the evolution by a time to of a 
point T for which 4> to j 0+T takes the same value of <p ^ T (X), i.e. 



0o,r(^°r) 



^ S s (S to T)) ds 



t +T 



to 



cb(S s T)ds = ^ t0+T (T) . 



Thus, Eq.(31) expresses the time evolving probability of the fixed set 4>o tT \( a ,t 

Mto(^0,rl(a,6)) = ^t {S to ^ t0tto+T \(a,b)) ■ 



(36) 



(37) 



Provided convergence to the steady state occurs, this equation yields a more and more accurate ap- 
proximation to the steady state measure of 0otI(o,6)i as *o grows, at fixed r. 

Equation (24), which is expressed in terms of the initially chosen phase space probability distribu- 
tion, and of evolving phase space sets, can now be given in terms of the evolving PDF, since 

Mto(0o,rU+) ^{S to 4> toM+T \ A +) Mo(0t O)to+ rU+) 



Mt o (0O,rU-) ^o( 5 *°^ ,to+rU-) W)(</Wo+rL-) 

This produces what we call 
(A-FR: 



(3f 



(exp(-^)>^ 1 



+ 



(39) 



Relation (39) is exact and holds for all to,r and possible pairs A and —A (whose range must be 
determined case by case). Letting (j> t0tto+T = £lt ,t +T and rearranging, Eq.(j27j yields the 

ft-FR: 

A + e(5, t , A,r) - - In ( e - n °.*o- n 



1 M*o(^0,tU+) 
-In = — 



t +r,2t + 



*0>'0 + 



(40) 



where the conditional ensemble average in the last term is computed with respect to the initial distri- 
bution /. 

Let us consider some direct consequences of the above relations. For A = 0, 6 > 0, to,r G M, and 
any odd </>, the conditional average of Eq.((39|) yields 



(exp (-Sh,t)h 



/^oC^o.tIc-^)) 

M*o(0O,t I (-*,«)) 



(41) 



which, in the 5 — ► oo limit, leads again to the NonEquilibrium Partition Identity (17), for the full 
ensemble average. 



Taking A = in Eq.(40), one obtains 



t ,t +Te(-<5,5) 



3 e(5,t ,A=0,r)r 



(42) 
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which shows that this particular conditional average is uniformly bounded in to, and tends to 1 when 
5 — ► 0, at fixed r. Equation (41 1 with <f> = O differs from (42 1 because of the integration times in the 
argument of the exponential function. For fixed 5 and r, such that St <C 1, one can write 



/g— ^0,t ~^t + T ,2t +T \ 



^t ,t +r6(-5,5) 



1 



(43) 



even for to - ► °°, which is interesting for the validity of the steady state Q-FR. Indeed, the last term 
in Eq.(40 1 ought not to dominate over the others, if the steady state f2-FR is to hold. Exchanging the 
roles of Af and Aj in Eq.(40), one obtains 



a e(S,t ,A,T)+e(S,t ,-A,T) 



(44) 



'*0' t 0+ rt ^ i (5 



r6A+ 



for any 5 > 0. 

One last property that is worthwhile to consider is the following asymptotic result, which is par- 
ticularly close to a steady state FR. As it happens in many situations of interest, let |0| be bounded 
by some 0* > 0, so that: 



In ( e - n o,t - n t +T,2t +T^ 



for all r. Then, taking 5 < 7, which implies |e| < 7, and 

T[t ) > = T (t ) , 



< 2t ^* 



one obtains the following 



7 



1 , Mt(^o,r(t)U+ 
A - 7 < hm — - In — = - 

t-oor(t) n t (n {t) \ A 



(45) 



(46) 



< A + 7 



(47) 



Recalling the connection between Q and the physical dissipation, which holds under proper choices of 
the initial /, the above condition can be written as 



r > 



2t J*F e V 



7 



(48) 



where J* is the maximum of the flux. This shows that, differently from the case of the steady state 
A-FR, the convergence time of the f^-FR does not increase when the dissipative field tends to 0; to 
the contrary, it may decrease. 

As in the transient case, an asymptotic 0-FR holds for A in a range which depends on the dynamics 
and on the observable, but for large r it remains close to (— £T , Q*). Indeed, one has 



2t + r 



1 /M0O,rUj 



ft* < - In 



< 



2t + r 



n* 



(49) 



and taking any 7 > and letting t grow, with r(i) 3> t, one can write 
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</>oo-FR: 

i M</v(t)L+) „ , N 

- O* - 7 < lim — In — _ 6 < ft* + 7 , (50) 
*-°or(t) Mt(0o,r(t)U-) 

for any odd 0, and with the smallest U*, among the upper bounds of the allowed dissipation functions. 
Because at equilibrium the dissipation function given by Eq.0) is O = -F e JV/k B T = 0, one has 

1 Mt (<?VL+) , N 

-7<~ln - .* <7, (51) 



whatever odd variable <\> is considered, SJ, A or J for instance. This yields the usual equilibrium 
symmetry, in the full range of J, not just in 0, and without the need for long convergence times. 

Given the dynamics S T , either there is a unique allowed density /o, hence a unique dissipation 
function Q or, given a second initial density fo, there is a corresponding O / Q. If this is the case, and 
the attractor for both fo and fo is the same, one can write 



o,t\a 



exp -n 0tt )„ » — \ (exp (-fto,*)}* ej4 + > ( 52 ) 



for large to, all r and all allowed A. Then, as far as the asymptotic fluctuation relations are concerned, 
the different dissipation functions are equivalent. 



Let us note that Eq.(47) does not necessarily imply the steady state O-FR, even when there is a 



single steady state and all initial measures //q with density fo are attracted by the same /j,oo, because 



the limit in Eq.(50) could be singular. Although this is not impossible, it is not the typical situation 



of particle systems of physical interest (cf. Section 4 and Appendix B), and Eq.(47) is a result rather 
close to the steady state J7-FR. 



In Section 4 we discuss the conditions under which Eq.(47| describes the statistics of trajectory 
segments sampled from a single steady state trajectory. Before we do this, we note that whereas in 
Sections 2 and 3, the initial distribution fo can have arbitrary form, provided it is ergodically consistent 
with the dynamics considered, a further condition on fo is imposed, in order to obtain the steady state 
J7-FR: the decay of the ^-autocorrelation, with respect to fo. This is necessary for the initial state 
to relax to a steady state. The further requirement that the ensemble is not only invariant under S , 
but that is also generated by a single trajectory, has the purpose to avoid the presence of distinct 
steady states, selected by different choices of the initial conditions. These states, indeed, would have 
different properties, and a unique fioo would not be possible. For instance, if fo was selected to be a 
uniform distribution (so that O = —A) and isokinetic dynamics was carried out, a distinct /Uqo would 
be expected for each distinct value of the kinetic energy. 

If the system under investigation does not converge to a steady state with positive and negative 
fluctuations, there is no need for a steady state relation, but the fi-FR, the 0-FR, the Oqo-FR and the 
0OO-FR remain valid for the given ensemble of trajectories. 

4 STEADY STATE FLUCTUATION RELATIONS 

We assume the steady state exists and that it is unique. This will not always be the case, for example 
if we start from a canonical phase space distribution at time zero and have isokinetic dynamics subse- 
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quently, there will be distinct steady state attractors for each trajectory with a different kinetic energy. 
Each of the steady state attractors will have different properties because they are at different kinetic 
temperatures. By saying that the steady state is unique we mean that, apart from a set of measure 
zero, the statistics generated by each trajectory is independent of the initial phase, in the relevant 
equilibrium phase space. 



Equations (27) and (40) are both exact, but contain information on the trajectories before the 



steady state is reached. To be of practical use in the description of the steady state fluctuations of the 
dissipation function, the last two terms in the equation have to be negligible. Therefore, for A ^ 0, 



some assumption on the nature of the logarithmic term of Eq.(40) is necessary since, in principle, this 



term could be of order O(to/r), and to should tend to infinity, in order to reach the steady state, before 
r does [2S]. We argue that, in most cases of interest, the dynamics of particle systems used to model 
nonequilibrium fluids satisfy the following: 

Steady state Jl-FR. For any tolerance 7 > 0, there exists <5 7 ,r 7 > such that 

1 Moo(^0,tL+) 

A- 7 <-ln s - <A + 1 (53) 

8 

holds for < 5 < <5 7 and r > r 7 , if both the values A and —A are (^-possible. 
Here, the following definition has been used: 

Definition. The value A E M is called <5-possible if given 5 > 0, there exists ta,& > such that 

^oo(^0,rU+) > 

6 

for all t > r A ,sf\ 

It follows that, for A to be 5-possible, /ioo(^o,r|^+) must decay with r more slowly than the error term 

e(5,t ,A,r) - \ In (e^^-^o+r^o+r^ 



of Eq.(40). Indeed, given 5 > 0, a pair (—^4,^4) 



and a fixed r, the logarithmic term in the right hand side of Eq.(40 ) may either take arbitrarily large 
values or be bounded, in the to ^ 00 limit. In the first case, because A, 5 and r are fixed real numbers 
and |e| < 5, the corresponding divergence is reflected on the measures ^j (Oo,r|^+)) and /Xf (Oo,TL-) 5 

8 S 

one of which at least either tends to zero or does not converge at all. Indeed, if for every R, r, to > 
there is to > to such that _ 

/^o(^0,tU+) 



R < 



one has 



R < 



/^o(^0,t| a -) 



or 



Mt ( n o,rL-) < 



1 



8 Because of the finiteness of r, this can be the case even though Hoo(£Io,t\ a +) 
state ensemble average of fi, (fi)oo say, does not belong to Af . 



R 

0, for r 



(54) 



(55) 



00, as it does if the steady 
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which corresponds either to ^too(^o,r 



0, or to the fioo non-measurability of Oo, T L-- In either 



case, —A is not J-possible, and similarly one may treat A. 

Clearly, A is of no interest for a steady state FR, i.e. it does not belong to the domain of the steady 
state O-FR with tolerance 7, if A and —A are not both <5-possible; in a number of cases that domain 
may even be emptyj^] Therefore, we focus only on the A values with bounded logarithmic term in 
Eq.(40), and let M(A,6,t) be the real number to which the conditional average in Eq.(40) tends, in 
the to — > 00 limit. We may then write 



A - - In M(A,S,t) 



1 A t oo(^0,rL+) 

6 < - In = ^- < A 

T jUoo(^0,TL-) 



1 



}nM(A,S,T) + S 



(56) 



If (1/t)\iiM(A,5, t) tends to zero, for the given A,S and growing r, the steady state Q-FR holds for 
A. This is the case if correlations decay, so that the thermodynamic behavior sets in. In particular, if 
the O-autocorrelation function is a <5-function (as e.g. in Ref.[MlEZ]), one can write: 



^g-^0,i -^t + T,2t +T^_ 



t 0' t + 



and 



+ = ^ e ~ n o,t -nt +T,2t +T^ — ^ e ~^o,t ( )^ ,/ e -f2t +T,2t +T^ 



which implies 



3 -n s ■ 



1 for all s, t 



(57) 

(58) 
(59) 



because of the NonEquilibrium Partition Identity (17). Then, (1/t) In M(A, 5, r) identically vanishes 



for all t. This is an idealized situation and, indeed, the tests performed on molecular dynamics 
systems, 10 cf. Ref . |4"U] . indicate that typically there exists a constant K (often not too far from 1), such 
that 

< — < ( e - Q <U0-^0+-.2i0+A < K ( 60 ) 

which may be understood as follows. In the first place, the conditional average is limited for growing 
to, if —A, A are ^-possible, and equals 1 for r = 0. Then, if the ^-autocorrelation decays with time 
scale tM^one expects InK to be of order O(ijvf)- Taking r sufficiently larger than tM, and letting 
ca be a natural number dependent on the given A, one may write 



< 



t +T,2t Q +T 



t 0> t + 
t -c A t M 

-Clo.tn-CAt 



-CI 



t Q +T + c A t M ,2t +T . 



n t0 ,t +T&Aj 



g - 2 tg+T + C a t M ,2*Q+T 



^t ,t +reA+ 



g --u,'-o~ c A t M 



(61) 



9 There are various reasons for which A or — A might not be in the the domain of the steady state f2-FR. For example, 
A could be larger than the largest val ue of fi, or there could be no fluctuations in the steady state (cf. A ppen dix B). 



The conditional average of Eq. ( 40 1 can, of course, be computed independently of the other terms in Eq. ( 40 1 , although 



this is not strictly necessary, because Eq.(40l is an exact identity. 

n This time is called Maxwell time; its order of magnitude is that of the mean free time. The Maxwell time expresses 
a physical property of the system and depends only mildly on the external field; usually, tM(F e ) = £m(0) + O(F^) 
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with an accuracy which improves with increasing to and r, if ca is the number of Maxwell times 
necessary to obtain the desired decorrelation at the chosen A. Indeed, one would eventually obtain 

to,T > CAtM- 

If Eq.(61) holds, the t —> oo limit of (e~ Qt o+ T+c A i M,2*o+^ i s finite for all <5-possible pairs (—A, A), 
and is insensitive to variations of r above a certain threshold. Indeed, let f = r + c^m and observe 
that 



H +f,2t + T 



f > 2t o+-( r )/ (r)dr 



(62) 

is the average of e Q °' t o- c A t M computed with the probability distribution obtained from the evolution 
of the initial one, for the time t = to + f. The value of this average may be estimated noting that, for S- 
possible pairs {—A, A), and growing to, it converges approximately to the positive number M(A, 5, f ), if 
Eq.|6l| holds. Furthermore, for fixed to > and for t of the order of the Maxwell time, (^e~ n °' t o' c A t M ^ 
approaches its steady state mean value, °'*o- c a*m ^ which is equal to 1 at equilibrium for all to 
and ca- Then, for sufficiently large to and t = 

M(A,S,t) ~ ( e - Q <M0-^0+r,2i +r\ 

(63) 

thanks also to the fact that f is by itself larger than the decorrelation time. Because the left hand 



to + t, one may write 



side of Eq.(63) does not depend on to while the right hand side does not depend on r, one concludes 



that further growths of to or of r do not cause any substantial changes in the conditional average 
in the middle. Close to equilibrium, this term should only have a correction of order O(F^) to its 
equilibrium average, 1. Numerical studies show that this is precisely what happens to the quantities 
in Eqs. Q60|61|62 |), 

Note that ca may vary with A, as the set of trajectories involved in the conditional average changes 
with A; ca could then even diverge. However, close to equilibrium at least, the existence of the transport 
coefficients implies that the conditional correlation functions must not grow too much (cf. Appendix 
B). 

Note also that, for the steady state Q-FR to hold, it suffices that (1/r) 1ilM(A, S, r) can be made 
as small as one wishes, by taking sufficiently large r and small 5. Thus, it is not necessary that 



Eq.(60) holds; it suffices that M(A,S,t) grows less than exponentially fast with r, or even that it 
grows exponentially fast, with a rate which tends to zero when S does, as in Eq.p2|), where ^4 = 0. 



However, the numerical tests performed so far indicate that Eq.(60 ) is usually satisfied in the full range 
of numerically or experimentally accessible A's. 



If the ^-autocorrelation decays, and Eq.(60 1 holds with the condition on O replaced by an equivalent 
one for another observable 6, one obtains the 



Steady State 0-FR. For sufficiently large t, 



lim 

to— »oo 



ln(e-^o^) 1 



1 Voo{<Po,t\a+) 

7 < — In = — < lim 



ln(e 



-a 



*o^o+ T 



-1 



+ 7 (64) 



This expression states that the ratio of the probabilities that (ft Q £ A$ and ^ 0iT G A s is determined 
by the average of the exponential of the values taken by O in the same time intervals in which cf) has 
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average in Ag. Because 4> can be any odd observable, O appears to be a rather special function of 
phase. 

One may wonder how the above could be reconciled with the modifications of the steady state 
A-FR suggested for a special kind of non-transitive dynamical systems (55J E5J As explained in 
Appendix B, Eq.(60) is not in contradiction to the scenarios of Refs.[2"5 j because the decay of 

correlations enjoyed by their dynamical systems concern the invariant measures, while the conditional 
average in Eq.(60) is computed with respect to the initial measure. Therefore, the arguments of 
this section and those of Refs.[3SJ EE] are compatible, since Eq.(60) would not hold if the cases of 
Refs.[2H[SIEII2I] were realized. The scenarios under which the approximation (60) does not hold, appear 
however quite peculiar, and this section explains why quite commonly the steady state O-FR does hold: 
time reversibility and the decay of the ^-autocorrelation characterize the most common deterministic 
models of nonequilibrium statistical mechanics. Through the connection with the Maxwell time tjv/, 
this approach also justifies why the convergence times of the steady state fi-FR do not grow with 
decreasing dissipation: the fact is that the O-FR depends only on material properties of the system, 
and these do not change much around equilibrium. This is quite different from the case of the A-FR, 
whose convergence times increase as the driving is decreased, and typically diverge in the equilibrium 
limit. As stated above, the calculations of this section hold for any regular density /, i.e. for all possible 
dissipation functions O, but the Q-FR turns into the S-FR only when the initial distribution can be 
generated by single field free trajectories, and AIT holds. 



5 CONCLUSIONS 

In this paper we have investigated the mechanisms underlying the validity of the transient and steady 
state FR's, following the prescriptions given by Evans and Searles over the last decade. As already 
evidenced in the literature, the transient relations hold under very general conditions; it suffices that the 
dynamics of the particle systems at hand be deterministic, autonomous, and time reversal invariant. If 
the initial state is properly selected, the corresponding Q, equals the energy dissipation rate S. Different 
initial densities / are allowed, for any given dynamics S T ; they are not completely equivalent, as in 
general they generate different dissipation functions $7, but to some extent they are interchangeable 



(cf. Eq.(52|). The transient relations describe the properties of evolving ensembles of systems, even in 
their asymptotic form. 

The steady state fluctuation relations need further assumptions. Firstly and quite obviously the 
system must relax to a nonequilibrium steady state. This requirement is nontrivial. The necessary and 
sufficient conditions for this to occur are unknown, but necessarily involve some loss of memory of the 
initial state, and the corresponding decay of correlations. Certainly there are examples of dissipative 
particle systems in contact with a thermal heat bath that do not relax to a steady state. Also, there may 
be distinct basins of attraction in M with distinct steady state measures. For example if M is taken 
to be all of phase space with an energy less than some bound and the initial ensemble is canonical, 
but the dynamics is isoenergetic, then obviously there will be no mixing between trajectories that 
have different initial energies, and the subensemble members will generate distinctly different energy 
dependent steady states. Different steady states may coexist also when dissipation is too high |H] . 
Transitivity (i.e. that a typical trajectory explores the whole phase space relevant to the given dynamics) 
avoids these difficulties. 
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As a matter of fact, most of the molecular dynamics studies performed so far indicate that conver- 
gence to a steady state independent of the initial phase is a very common situation. This convergence 
is consistent with laboratory experiments where for systems with low Reynolds/Raleigh numbers the 
nonequilibrium steady state (if such a state exists) is quite independent of the initial starting state or 
the preparative history of the system. The low shear rate viscosity of fluids is empirically observed 
to be a function of the temperature and density (or pressure) of the fluid. The independence of the 
initial phase or preparative history is the nonequilibrium analogue of the situation where at equilibrium 
thermodynamic free energies are commonly observed to be thermodynamic state functions. Even at 
equilibrium the conditions required for this to be the case are not known (they are violated in glassy 
systems for example). 

If a single steady state is approached, its statistics may be recovered from an individual trajectory, 
and this may have no knowledge of the initial ensemble /. This loss of knowledge of the initial ensemble 
amounts to a decay of correlations, which may be further weakened. It must be stressed that the decay 
of correlations required by the validity of the O-FR reduces to the ^-autocorrelation decay, with respect 
to the initial ensemble /. As is well known, this is required, close to equilibrium, for the existence of 
the transport coefficients. 

In the case that transitivity and boundedness of f2 suffice to obtain the steady state Q-FR (some- 
thing that needs further study), we observe that these are much weaker conditions than the Anosov 



assumption, which also implies the boundedness of the phase space contraction rate _ This is in- 
teresting, because one can name a number of cases which do have bounded like the isokinetic 
(or isoenergetic) hard particle systems; the isokinetic or Nose-Hoover thermostatted systems for heat 
current and shear flow systems with bounded interaction potentials; the color conductivity isokinetic 
systems with unbounded potentials etc. Transitivity, in turn, is almost always verified. 

We note that in order to obtain a steady state S-FR, the initial distribution, /o, must be obtainable 
from a single the field free evolution, and AIT must hold. These conditions are easily obtainable for 
systems of physical interest. 

The E-FR's directly refer to the dissipative fluxes or to quantities simply related to them, rather 
than to the phase space contraction rate. This avoids the difficulty of justifying how a A-FR can speak 
of the fluctuations of the entropy production rate. As discussed in [5], these difficulties are more evident 
close to equilibrium, where one observes that the convergence times of the A-FR typically diverge as 
the dissipative applied field decreases and the system approaches equilibrium. This is due to the fact 
that A contains undesirable terms, which have to average to zero. As one approaches equilibrium the 
dissipation decreases as the square of the dissipative field while the undesirable terms remain constant 
(to leading order in the field). Therefore, longer and longer averaging times are required before the 
A-FR converges. It has been argued [24J that the S-FR can be obtained from the A-FR because the 
time averages of £ and —A become equivalent in the long time limit. However, this argument can 
not explain the satisfaction of the S-FR at physically accessible times in non-isoenergetic systems at 
low fields. Moreover, if the Van Zon - Cohen picture applies to deterministic particle systems, 
as argued in Refs.0 [371 32], the asymptotic equivalence between S-FR and A-FR is restricted to a 
limited fraction of the possible fluctuation values. As a matter of fact, E equals —A only if / is uniform 
in M, which requires M to be compact. When this is not the case, £ and A are substantially different 



12 Recall that the Anosov assumption was considered appropriate for the particle system of |B], in which —A = fl, 
although was not bounded. 
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quantities, as even Ref.[3Tj shows. The difficulties with the weak field convergence of A-FR make the 
derivation of Green-Kubo relations from A-FR, highly problematic [5]. Differently, if a relation for 
the dissipation is obtained independently, this difficulty is not present. In particular, the convergence 
times are understood in terms of material properties related to known, physical properties of the system 
(e.g. the Maxwell time is related to the infinite frequency shear modulus and the zero shear rate shear 
viscosity), and the derivation of the Green-Kubo relations from the steady state S-FRs is reasonably 
straightforward [S]. 

It is also important to consider the consistency of the results obtained in this manuscript with 
those obtained for the axiom C systems postulated in Refs.[55| 155]. As discussed in Appendix B, these 
systems may satisfy the decay of correlations of all observables, with respect to an invariant measure 
supported on a surface M r , of dimension lower than that of the phase space M. Then, assuming as 
in Ref.[M], that A is proportional to the contraction rate of volumes in M r , A r say, A should obey a 
modified FR. So far there is no direct evidence, in particle systems, of the validity of that assumption 
(supported by numerical studies of certain hydrodynamic models [21]). In fact, in the case considered 
in Ref . [T2] . the dissipation is sufficiently high that the dynamics is not transitive yet the S-FR holds 
without the modifications conjectured in [32]. In other words, our work and Refs.JSSJ [36] usually 
consider different observables, and when they consider the same observable, it is not necessary that 
they lead to the same conclusions of Refs.[5S| 136]. If, on the other hand, the treatment of Refs.[3SJ |36J 
applies, the decay of correlations of Section 4 is compatible with it, as explained in Appendix B, 



because it is referred to the initial measure, and not the invariant one _ The energy dissipation E 
remains related to A, not to A r , even if the steady state is confined inside a lower dimensional manifold 
M r . 

As pointed out by one of the anonymous referees, the existence of a unique stationary measure 
may be considered a strong assumption, and the necessary and sufficient condition required for it to 
hold are not known. Therefore, the study of Anosov systems constitutes a natural starting point for 
the purpose of classifying the dynamical systems that verify some kind of FR. The purpose of the 
present paper is different: it amounts to analyzing systems which cannot be chosen at will, but are 
selected by physical conditions, in order to understand further their properties as well as the physical 
mechanism responsible for FR's to hold in natural systems. We have found that time reversibility and 
ergodic consistency imply, without further assumptions, a number of transient and asymptotic FRs 
(cf. Sections 2 and 3), which are amenable to experimental verification. If the system of interest does 
not reach a unique steady state, assuming that it is Anosov does not help of course, but those relations 
remain valid. When present, convergence to a unique steady state is the manifestation of a certain 
decay of correlations of the observables, which typically are required to relax to a steady state. The 
approach proposed in this paper thus makes it clear that the details of the microscopic dynamics (like 
the very strong hypothesis of uniform hyperbolicity) do not play any role in the validity of the FR's, as 



appropriate for thermodynamic relations. The class of dynamical systems which obey Eq.(60) cannot 
be determined, however, by this approach. 

These facts may help in developing further the theory of nonequilibrium phenomena, as the deriva- 
tion of the novel relations given above indicates. 



13 One should also observe that Eq.(40i is exact, hence cannot lead to any contradiction about reversible dynamical 
systems. 
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APPENDIX A 



Using the general definition of the dissipation function (13), it can be shown that the instantaneous 
value of the dissipation function is directly related to the instantaneous dissipative flux for a wide 
class of systems. Here we consider three cases: isoenergetic, isokinetic and Nose- Hoover thermostatted 
dynamics. 

A general form of the equations of motion for a field-driven iV-particle, thermostatted nonequilib- 
rium system is 

q; = Pi /m + Ci(r) F e 

Pi = Fi + D;(r) • F e - a Pi (65) 

where F e is the applied field that is coupled to the system via the phase functions C;(r) and Dj(T), 
The term — api is a deterministic, time reversible term used to add or remove heat from the system 
The adiabatic equations of motion 

qi = pJm + Ci(T) F e 

Pi = F 4 + D,(r)-F e (66) 



simply lack the thermostatting term. If Eqs.(66) are Hamiltonian, the phase space expansion rate A 
obviously vanishes, a condition referred to as adiabatic incompressibility of phase space, AIT [29J, but 
AIT may hold even if the general adiabatic equations are not Hamiltonian. Here, for simplicity, we 
concentrate on systems which satisfy AIT, because their dissipative flux and dissipation function are 
simply related. In the other cases, the relationship between dissipative flux and dissipation function 
is only more complicated. We also assume that the initial distribution of phases fo may in principle 
be generated by a single field-free (F e = 0) thermostatted dynamics. By this we mean that fo must 
not only be preserved by the equilibrium dynamics, but also that a single field free trajectory may 
explore all of its support. In other words, that this support is not the union of disjoint invariant sets, 
like the constant energy surfaces corresponding to different initial energies of Hamiltonian systems. 
This is consistent with nonlinear response theory that shows that if one is interested in the statistics 
of the physically relevant steady state observables, fo should be that which corresponds to the field 
free dynamics (2HJH3]. The possible difficulties that might be encountered if fo is not selected to be 
the distribution function generated by the field free dynamics can easily be seen by considering the 
distribution of 4> Q T along an equilibrium trajectory. In this case the correct physical picture requires 



the right hand side of Eq.(39) to be 1 for all A, to, 6 and r. This can only be the case if fl(T) = 
for all T. Using Eq.(34), this implies that fo(SQT) = f^S^T) under the field free dynamics Sq, and 
that fo = ft must be preserved by the field free dynamics. Note that a single trajectory can't possibly 
generate the ensemble in physically relevant times, in any system made of more than a few particles, 
but it may do it, in principle, allowing it to evolve for exceedingly long. This is needed for the possible 
initial equilibrium state to be uniquely determined. 

In general, the dissipative flux, J, is obtained from the adiabatic time-derivative of the internal 
energy, Hq = YliLi + ^Kq)) 1 - e - from the derivative of the internal energy under the dynamics of 



Eqs.(66): 



N 

H$ d = ~ ■ D 4 - Ft ■ Cj) • F e := -JV ■ F e , (67) 
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while the full dynamics of Eqs.(65) yields 



N 

H = -3V-F e -aJ2 



Pi -Pi 

m 



-JV ■ F P - 2Ka 



(68) 



were K is the kinetic energy of the system. Consider now three widely used deterministic, reversible 
thermostatting methods. 



Gaussian isoenergetic (ergostatted) system. In this case Hq is fixed, i.e. Hq = 0, and Eq.(68) 
gives 

J(r)y-F e 

° (r) = " 2K{T) {69) 
For the ergostatted system described by equations ( |65| ), the phase space expansion rate is, 

A(r) = V • f = ^ ( dC J^ + g "l (r) l • F e - dNa(T) + O n (1) A = -dNa{T) + O n (1) (70) 



i=l 



where AIT is used to obtain the final equality, and Ojv(l) is a correction which is order 1 in N. 
The Otv(1) term can be explicitly determined from the partial derivative of the Gaussian ergostat- 
ting term with respect to p. The equilibrium phase space distribution function is microcanonical 



on a surface Hq =const, hence f(iS T T) = f(S T T) = f(T), and Eq.(13) yields: 



njf (r) = -A 0)T (r) 

for all r, so 

n(T) = -A(r) = dNa(T) + O n (1) = - J(T) • F e [3(T)V + N (1) 
where 0(T) = (dN-d- l)/(2AT(r)). 



(71) 



(72) 



Gaussian isokinetic system. In this case the kinetic energy is fixed K(T) = Kq = (dN — d 



l)ksT = (dN — d — 1)/P- Equation (68 1 can therefore be used to show 



a(T) 



(H (T)+J(r)V-F e )(3 
dN-d-1 



(73) 



while the phase space expansion rate A is still expressed by (70 ) to order N, with a different Ojv(1) 
term. The equilibrium phase space distribution function takes the form f(T) ~ e~^ H °5(K(T) — 
K ), on a surface K(T) =const, hence /(r)//(i,S T r) = e^fo «o(r(a))<fa ) and Eq.O leads to: 



O 0)T (r) = p / H (T(s))ds - Ao, T (r) = -(J • F e ) r VP + Ojv(1) 

JO 



(74) 



for all r. Combining this with Eqs.(73) and (170]), one obtains 



E(r) = -j(r) • F e vp + Ojv(i) 



(75) 
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• Nose-Hoover thermostat. In this case, equations (65) are supplemented with the equation, 



a = ^{2K{T) - dNk B T) (76) 

where Q is a constant related to the relaxation time of the thermostat, and T is the imposed 
average temperature. In this case, 

a(T) = -GHi,(r) + J(r)V • F e )/(2K(T)) (77) 

and 

A(r) = V • r + ^ A = -dNa . (78) 
oa 

The equilibrium phase space distribution function is the Nose-Hoover extended canonical distri- 
bution f(T) ~ e -/3(^o + |Qa 2 ) ; hence /(r)//(z5 T r) = e Pfo(H T(s)+Qaa)ds_ g Q equation (g) can 

be used to show that: 

(r) = (3 f T (H (s) + Qaa)ds - A 0>T (r) = -(J^F^ T V7? (79) 
■> o 

for all r, and then 

E(r) = - J(r) • F e V(3 (80) 

Therefore, for field driven nonequilibrium systems whose equations of motion satisfy AIT, the dissi- 
pation function and the dissipative flux are simply related. In fact, for thermostatted systems with a 
constant field, S(r) = kJ(T) where fc is a constant. 



It is straightforward to show that the same expressions (72), (75), (80) can be obtained even if only 
a fraction of the system's particles is subjected to a thermostat (e.g. the wall particles), see [11J,|12J 
for example. Furthermore, if the walls are large, the results are not sensitive to the details of the 
thermostatting mechanism |13j . 

APPENDIX B 

The derivation of the steady state relations from transient ones needs a careful analysis of the condi- 
tional average 

( e -nO,*0 • e -^0+r,2t +r\ (81) 



which appears in Eqs.( |27|40~ ), and of its limit 



M(A,5,t)= lim (e^ ^ ■ e~ nt o +T * 2t o+ T ) rj c .+ (82) 



to— »oo 



which is needed for the transient states fit to eventually reach the steady state /Xqo. In particular, 
M(A,5,t) must exist and be positive for r larger than a certain ta,s, for the pair (—^4,^4) to be im- 
possible, i.e. for fluctuations of size close to A to occur, and for A to be in the domain of the fluctuation 
relation with some tolerance 7 > 5. 
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As observed in the text, Eqs.(27) and (40) are exact, hence cannot contradict any correct result on 
deterministic, reversible dynamical systems. The question is whether one can extract sufficient infor- 
mation from them, that they can be of practical use for the description of the statistics of steady state 
trajectories. This requires that either the dynamics of interest be explicitly given, or that the possible 
situations be explored. In section 4, we have considered the possibility that the ^-autocorrelation 
decays in time, as the most common for the deterministic reversible particle models of nonequilibrium 
fluids. In fact, less is needed for the steady state O-FR to hold for A; it suffices that M(A, 5, r) does 
not grow too fast with r: subexponential growths, and exponential growths with rate not larger than 
5 (as for A = 0) are all acceptable. 

In the linear regime, the decay of correlations in the steady state can be well approximated by 
their decay at equilibrium, and for this reason, in the linear regime, the transport coefficients are 
given by the integral of the equilibrium time correlation function, multiplied by the field |29j . Lack 
of decay of correlations of f2 corresponds to the non-existence of the transport coefficient associated 
with Q. Further away from equilibrium, nonlinear response theory applies. Now the time correlation 
function is still determined with respect to the initial measure fi, but the dynamics used to compute the 
nonlinear transport coefficients is the nonequilibrium one (cf. the Transient Time Correlation Function 
|29J). This is analogous to what is done in calculations of M(A,6,r), and if the nonlinear transport 
coefficients exist in the nonlinear regime, then there should be a decay of the ^-autocorrelation, with 
respect to the equilibrium measure. Note also that the growth of r increases the separation of the 



end of the integral in the first exponential of Eq.(81) from the beginning of the integral in the second 



exponential, hence that the growth of r should contribute to the decorrelation. 

Therefore, it looks quite plausible that, in most cases, one of the above considerations on the decay 
of correlations be satisfied. There is, however, one competing effect, which may become important in 



particular circumstances: the set of trajectories over which the conditional average (81) is computed 
changes with r, and this may balance the decorrelation phenomena in some cases. 

To illustrate these facts, consider a simple model that has been discussed in connection to the 
fluctuation theorems in the past [2Ej; and apply the procedure described in this manuscript to that 
example. The model consists of a particle moving in empty space, under the action of a constant 
external force F e and a Gaussian thermostat: 

q = p 

P = F e p 

p p 

where the kinetic energy K = p • p/2 is a constant of motion. The dynamics of this particle is rather 
simple: for initial conditions with p pointing in the direction opposite to F e , p is a constant of motion; 
for all other initial conditions the direction of p tends to the direction of F e , while the magnitude of p, 
V2K, is constant. The treatment of Sections 2 and 3 applies to this system, since it is reversible and 
an initial distribution that is ergodically consistent with the final distribution can be selected. However 
there are various reasons why the treatment in Section 4 cannot be applied to produce a steady state 
S-FR in this case. For instance, the distribution generated by the field free equations of motion is not 
ergodically consistent with the steady state dynamics (see Appendix A). In addition, the attractor is 
a fixed point in p-space, and the time-dependence of q is irrelevant; so, for any f2 that is a function of 
p only, there are no steady state fluctuations - that is for small S, no pair (—A, A) is 5-possible, hence 
that the domain of the steady state O-FR is empty. The theory outlined in this paper also anticipates 
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these facts, since correlations for any phase function do not decay, and therefore the term M diverges in 
the to — ► oo limit. In particular, selecting the dissipation function that would correspond to a uniform 
initial distribution, we obtain Q = —A = {d — l)F e • p/2K where d is the number of dimensions in 
configuration space. The range of values of this VL is [— (d — 1)\/2K |F e |, (d — l)y/2K\F e \] and any 
smooth probability distribution in it (e.g. the uniform one) evolves in time in such a way that 



O 



'0 if (d-l)V2K\F e \(£(A-5,A + 5) 

_ (83) 
k l if (d-l)V2K\F e \e(A-6,A + S) 



As all trajectories, except those with initial condition p(0) = — (d — l)v2^F e /|F e |, tend to have 
f2 = (d — l)v / 2KF e /|F e |, this means that the integrals Oo,t an d ^t +T,2t Q +r diverge linearly with to- 



Therefore, Eq.(60| does not apply in this case. This will be the case for any steady state that 
has no fluctuations. Nevertheless, the transient and asymptotic relations remain valid for the evolving 
ensembles. The transient Q-FR expresses the rate at which the initial probability of observing —A 
vanishes, compared to the initial probability of A, and the f^-FR expresses the corresponding asymp- 
totic rate; but the domain of the steady state O-FR is empty. A similar situation, with no fluctuations 
of opposite sign in the steady state, may be produced in more realistic models by applying extremely 
high fields, making the motion ordered. In these cases, the domain of the steady state fi-FR would also 
be empty. However, intermediate situations are possible, like those of Ref . [Ti] . in which the attractors 
have small dimension, but the initial correlations decay, and the steady state O-FR holds. 

References [SS1ES] consider systems with steady states confined within manifolds M r of dimension 
smaller than M., whose volumes fluctuate in time, so that the corresponding phase space contraction 
rate, A r , obeys the standard fluctuation relation. In these systems, correlations of all observables 
decay with respect to the invariant measure fioo, supported on M. r . Furthermore, if the phase space 
contraction rate A, in the space of the field-free dynamics M, is proportional to A r , as assumed in 
[55] . a modified FR holds for it, which can be expressed as 

1 Moo(Ao,tL + ) 

cA - 7 < - In = s — < cA + 7 (84) 

T ^oo(A ,r| A -) 

O 

with c < 1. The relevance of such a possibility for particle systems is not obvious, but no difficulty 
emerges with the theory developed here, since the decay of correlations with respect to does not 
imply the decay of correlations with respect to a non-singular measure supported on M.. Many different 
scenarios are possible, and should be analyzed case by case, but this goes beyond the purpose of the 
present paper. Here it suffices to observe that no scenario that can be realized can contradict the 
framework of the present paper; to the contrary, it could be explained in it, through the analysis of 



the conditional average (81). 
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